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This is a supplementary material of our recent paper [l], where a class of exactly solvable S = 
1 quantum Ising spin models were studied based on the hole decomposition scheme. Here we 
provide some details for the Green functions, the spin-spin correlation functions, as well as the spin 
susceptibility in the presence of dimerization. 



I. INTRODUCTION 

In Ref.|l[ we have studied a class of the 5 = 1 spin 
chains with the nearest neighbor Ising coupling and both 
transverse and longitude single-ion anisotropy by a com- 
binational use of a hole decomposition scheme and a re- 
cursive method. These models include the first example 
of the dimerized 5=1 quantum spin chain where all the 
eigen states can be solved exactly. In this supplemen- 
tary material we present some detailed derivations for the 
physical quantities of the 5 = 1 dimerized chain. All the 
notations are the same as in Ref.0. In Sec. EH we discuss 
the Green functions of the uniform or dimerized chains, 
respectively. In Sec. IIII1 we study the longitudinal spin- 
spin correlation function at zero- or finite-temperatures. 
In Sec. IVl and IVdl we list some detailed formulae for the 
segmented M-matrices and the partition functions. 



For the uniform system, the wavefunctions <i>kj and 
^kj can be taken as real, we have 



The Green function can be then expressed as 

GM =-$^tt*i*k,tanh|/9A(*)/2] . 

k 

Note that tanh[/?A(fc) /2] — ► 1 at the ground state (fi 
oo), so we have 

G jq (p^oo) = -J2*kj$k q ■ 
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II. GREEN FUNCTIONS 

A. Green functions of the uniform spin segments 

The original 5 = 1 quantum Ising model is mapped 
onto a large family of the segmented 5=1/2 transverse 
Ising models classified by the total number of holcs[2, 
3]. These segmented 5 = 1/2 models are then solved 
by introducing the Bogoliubov fermionic quasi-particle 
operators r)l and rjk as defined in Eq. (14) in Ref. [l[. 
Inversely, we have 



4 = E 



-% 



The Green function, or the two-point correlation func- 
tion, is defined by 



G jq =(Ff-)FM), 



(1) 



where Fj = cj ± Cj . 



We denote the wavefunctions for the chain with pe- 
riodic boundary condition (cyclic) and open boundary 
condition (free ends) by (<E> C , \l/ c ) and ($* , V I / ^), respec- 
tively. Then we have 



2/Zsinjfc , k > 0, 



M = 
kl A(jfe) 



2// cos jk , k < 0, 
D 



(3) 



[(1 + A cos k)$ c kj + A sin k$ c _ kj ] 



where I is the length of the segment. The Green function 
is 

G% = L r + XL r+1 , (4) 
where r = \ j — q\ and L r was defined in Refs. [3, Q 



L r = — / dk 



1 



Similarly, 



Vl + A 2 + 2Acosfc 



A k sin(j -q + l)k, 



: cos fcr 



(5) 
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where 



Ah = - 



21 + 1- 



sin (21 + l)k 



sin fc 



-1/2 



Consequently, we have 

G f jq = - A l 6 k sin i k sin 0' -i + l ) k - 



(6) 



(7) 



At the finite temperatures, we need to add the factor 
tanh[/3A(fc)/2] to Eqs.Q and 0. 



B. Green functions of the dimerized segments 

In the presence of dimerization, the wavefunctions <&kj 
and are complex in general. So we now have, 



F 3 ^X)*^* - *^*' 



(8) 



Then, the Green function is expressed by 

Where, (r^%) = [exp (A k /(k B T)) + l]" 1 , satisfying 
Fermi-Dirac statistics. At the zero temperature, the 
Green function can be written as 



G jq =D j Y\j,q}+2J j Y\j + l,q], 



(10) 



where 



,i(j-q)k 



-[l + (-l) J '+V7 



+ (-i)^ 7 + (-i)V]. (ii) 

The dimerization parameter 7 is defined by 

1 -r 



7 = 



1 + T 



(12) 



with t being determined by Eqs. (19) in Ref. [l[. 

Generally, r has two solutions, corresponding to the 
upper/lower signs of ± respectively in Eqs. (19) in 
Ref. [l[ . In order to numerically calculate the Green func- 
tion, we need to express Y[j, g]-function in terms of real 
variables. We introduce px,2, 91,2 to express complex 7 
as follows. 



7i=Pi+%j 72 = P2 + iq2, 



(13) 



px ; 2 and q\ } 2 are the real and imaginary parts of 71,2, 
respectively, 



Pl,2 
91,2 



b\ + &2 + 4M>2 cos2fc-(CiTC2) 
[(61 + 6 2 ) cos k - Ci ± C2] 2 + (62 - 61) 2 sin 2 fc ' 

-2(6 2 - 61) sin fc [(pi + b 2 ) cosfc + Ci T C2] 
[(61 + 62) cos fc - Ci ± C2] 2 + (62 - 61) 2 sin 2 fc ' 



where the subscript 1 corresponds to the upper case, the 
subscript 2 corresponds to the lower case. Ci 2 are given 
by 

Ci = («2 - ai)/2, 



C 2 = rVi -P + Qcos2fc . 

where ax, 02, Pj Q and T are defined in Ref.[l|. 

For convenience, we divide fc-region [— tt,tt) into two 
subregions: (I) for [— 7r/2,7r/2) and (17) for [— tt, -tt/2)U 
[7t/2, 7r), respectively. Thus Gj q can be expressed by 



— °jg "r" °j'<j ' 



(14) 



In Region (/), because of the symmetry between fc and 
-fc, the Green function can be reduced in (0, 7r/2), 



^ A_i(fc) 

(0,7r/2) U ^ 



{I>i[l + (-l) i+ *(p? + «?) 



+ (~l) J Pi + COs(j - q)k 

+ 2J 3 [l + (-iy + <i +1 ( P \+ql) 

+ (-l) j+l P i + (-l) q Pi] cos(j - q + l)fc} . 



(15) 



A similar Green function can be obtained for Region (II). 
The function Y[j, q] can be rewritten as 



Y\j,q] = 



^ A_i(fc) 



[i + (-iy^(pl + q D 



+ (-iy P i + (~l) q Pi] cos(j - q)k 



(16) 



(ir/2,7r) V y 



+ (-iyP2 + (-1)>2] COs(j - q)k. 

So it is convenient to express the total Green function 
Eq. (flU|) in terms of Y[j,q\. In the dimerization case, 
there are four such Green functions associated with the 
four different parity combinations of the segments. 



III. CORRELATION FUNCTIONS 
A. Zero temperature 

In this subsection, we discuss the spin-spin correlations 
at zero temperature. In Ref. [l[ we show that the ground 
state has no hole if D z > — A^(0), otherwise, it has holes 
once D z < — A^(0). In the latter case, the holes break 
the original chain into segments. We note that only the 
intra-segment spin-spin correlations are non-zero. 

For D z > — Afc(0), the spin-spin correlation function 
of S z is defined by C^ n = (* |5'm' S 'nl^ , o), where |* ) is 
the normalized ground state of the Hamiltonian. By use 
of the Jordan- Wigner transformation, one has 

C z mn = (VolFtfFMF^ ■ ■ ■ F^fM |*o). (17) 
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It is straightforward to show that (* |Fj ±) i ;l 9 (±) |* ) = 
iSj q . By further utilizing the Wick Theorem, we find 
that 



C z = 

^ 'ran 



G m +l,m+l G m +l,rn+2 



G n -i 



m+l 



G n -i 



m+2 



Gra,n 
Gm+lji 

Gn — l.n 



(18) 



for n > m, where, Gj q 



(* \F^F^\V ) 



(* |^ (+) ^ (_) |*o) 



The general expression of Gj q is derived in Sec. Ill Al 
for the uniform chain and in Sec. Ill Bl for the dimerized 
chain respectively. In general, one has 



G jq = D j Y\j,q] + 2J j Y\j + l,ql 



(19) 



where Y\j, q] is given by Eq. p6[) . For a uniform system, 
Y\j,q]=Y[q,j] = ±L j _ q . 



B. Finite temperatures 

At finite temperatures, the contribution from p ^= 0- 
sector should be taken into account. A recursion formula 
similar to Eq. (36) in Ref. [l| can be derived for the cor- 
relation function as following 



1 



Z(L) 



L L~p I 



p—0 1—0 m,n 



(20) 



Where, p z mn is the correlation function of individual seg- 
ments. It has a similar form with that in Eq. (fT5|) . but 
now Gj q should be replaced by Gj q ((3). 




FIG. 1: Temperature dependence of the spin-spin correlation 
function in a uniform spin chain with A = 1.5. 

In Fig. [TJ we plotted the temperature depen- 
dence of the spin-spin correlation function per site, 



E™,„Cw We find that when D z < -A h (0), 
the correlation function approaches to zero in the limit 
T — > 0. This indicates that the ground state is in 
the hole condensation phase. On the other hand, when 
D z > — Ah(0), the correlation function approaches to a 
finite value (about 0.85 for the two cases shown in the 
figure) in the zero temperature limit. 



IV. SPIN SUSCEPTIBILITY 

The spin susceptibility of the 5=1 QIM can be also 
calculated using the recursion formula introduced in the 
previous section. To do this, one needs to first evaluate 
the partition functions of each 5 = 1/2 Ising segments 
in the applied magnetic field £, denoted by z (£„,£). The 
partition function of the original 5 = 1 QIM is then given 
by Z(L,0 = £p=oE {M I&(^K- In terms of 
the segment magnetization m{l n ,T) = — -g ain ^"^) an(1 

the to- 



the segment susceptibility x{ln,T) 



dm(l n , T) 



tal susceptibility x(T) at zero-magnetic field can be ex- 
pressed as 

V ' p=0 1=0 



X {l,T)z{l)Z^- l \L-p-l). 



(21) 



Thus the hole decomposition scheme provide an al- 
ternative approach to calculate the susceptibility of the 
5 = 1 QIM. This approach is efficient provided that the 
susceptibilities of the corresponding 5 = 1/2 TIM's with 
varying chain length L are available. We note that the 
susceptibility of the 5 = 1/2 TIM has already been stud- 
ied by a number of groups 0, 0, HI]. So in principle these 
results could be used in the numerical study of the sus- 
ceptibility of the 5 = 1 QIM. 



V. DIAGONALIZATION OF THE M-MATRIX 

For a periodic spin chain, the diagonalization of the M- 
matrix has been discussed in Sec. IV A in Ref. [l|. Here 
we consider the diagonalization of this I x I M-matrix for 
an open spin chain with the length I. The aim here is to 
solve the following eigen equation 



M$ fe = A 2 (fc)$, 



(22) 



in various cases, where &k(j)'s take the form of Eqs. (23) 
in Ref. Q. 

We assume that the two ends of the open chain are 
located at the sites r% and r2, respectively. r% and 
can be either odd or even, so there are four kinds of M- 
matrices. In the following, we will present the results for 
each cases. 



4 



A. (ri,r 2 ) = 


(odd, even) 


The secular equation is 


In this case, the matrix M 


is defined by 


[{a 2 - a x ) ± W] [h sin(l + l)k + b 2 sin(l 


M = 


/a bi 
bi a 2 b 2 
b 2 a\ 


• • • \ 

• • • 

• • • 


26i(6i +b 2 + 26i6 2 cos2fc) sm(i + l)fc 

(23) 







■■■ ai bi 


where, 03 = D 2 . 




V 


■■■ bi a 2 J 





(29) 



where <xi,2) 6^2 are defined in the main text and clq = D\. 

The energy spectra can be solved following the ap- 
proach introduced in Section IV. The result is given by 



A2 ( fc ) - ^Xe-^ 1 ^- 10 ^ 

1 „ ( „2ik 1 — 2ik\ , l ^( Jiik 

+a 2 (e — t e e ) + b 2 T[e 
The reflection parameters arc 

f _ 2i(l+l)k 

t (bie lk 



t a e 



— 3ik\ 



(24) 



tc 



o 2 e 



{b x e- lk + b 2 e lk ) 

Then, the secular equation is given by 

[{a 2 - ai) ± W] [61 sin(; + 2)k + 6 2 sin/fc] 
2(a — »i)(6i + + 26i6 2 cos2fc) sin/fc 
= b 2 ' 



(25) 



where W is defined as in Eq. (20) in Ref. 0. 

Other cases can be solved by the same way and the 
results are listed below. 



B. (^1,^2) = (odd, odd) 
The reflection parameters t o e are 



f _ 2i(l+l)k 

t e (he- lk +b 2 e ik ) 
(bie ik + b 2 e~ ik ) 



(26) 



The secular equation is 

[(ai - a 2 ) ± W] [61 sin(Z - l)k + 6 2 sin(Z + l)fe] 
26 2 (6? + b\ + 26i6 2 cos 2k) sm(l + l)k 



a - ai 



(27) 



D. (ri,r 2 ) = (even, odd) 
The reflection parameters t o e are 

, _ 2i(i+l)fe 

t (b ie - lk + b 2 e lk ) 



(30) 



(6ie ife + 6 2 e- lfc ) 



The secular equation is 



[(ai - o 2 ) ± W] 

26i[6i sin(Zfc) + 6 2 sin(^ + 2)k] 



(31) 



a-3 - a 2 



VI. THE PARTITION FUNCTIONS OF 
SEGMENTS 



The partition function of individual segment of length 
I and parity (ri,r2) (defined in Sec.|V| is given by 



Z (r 1 ,r 2 ){l) = J| C0Sh 



feiG(0,7i-/2), 
fc 2 e(7r/2,7r) 



[j8Ai(*i)" 


cosh 


[/3A 2 (fc 2 )l 


2 




2 



. (32) 

where, k\^ 2 satisfy the corresponding secular equations. 



(ri,r 2 ) = ( 



even, even) 



The reflection parameters t a ^ e are 



2i(Z+l)fc 



t e (bie 



ik 



■b 2 e- ik ) 



(6ie" lfc + 6 2 e lfe ) 



(28) 
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